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Abstract. We carry out the microscopic calculations of the atomic EDM in 
deuterium, Xe and Hg atomic systems. Due to Schiff's theorem, we obtain the atomic 
EDM only from the finite nuclear size effects, but it is found that the atomic EDM 
is appreciably larger than expected. This is essentially due to the new mechanism in 
which the nuclear excitations are taken into account while atomic states stay in the 
ground state. The EDM of deuterium is found to be do — 0.017d„ while the EDM of 
Xe and Hg become dxc — 1.6d„ and dn g — — 2.8d n , respectively. It turns out that the 
new constraint on the neutron EDM becomes d n ~ (0.37± 0.17± 0.14) x 1CT 
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1. Introduction 

The existence of T-violating interaction in nature is one of the most fundamental 
subjects in field theory. There must be some T-violating interaction in nature as long 
as we believe the big bang cosmology since the big bang must have started from the 
T-violating interaction. 

In order to see the violation of the T-invariance, one measures electric dipole 
moments (EDM) of particles, nuclei and atoms in their ground states. Until now, the 
upper limit of the neutron EDM d n is around [H 121 IS] 

d n ~ (1.9 ± 5.4) x KT 26 e • cm. (1.1a) 

d n ~ (2.6 ± 4.0 ± 1.6) x 1(T 26 e • cm. (1.16) 

There have been many experimental efforts to measure the EDM of the atomic systems. 
The good examples are found for the EDM of 129 Xe [HE] and 199 Hg [7j. In this case, 
however, one has to be careful for extracting the EDM of electrons or nucleons from the 
measurement of the atomic EDM since there is Schiff's theorem |Hj. This theorem states 
that the EDM of the atom is canceled out due to the symmetry restoration mechanism 
as long as the constituents are interacting through the electromagnetic interactions with 
the nonrelativistic kinematics. In order to obtain the EDM of atomic system, one has 
to find the relativistic effects in the EDM operators Ej- Since electrons become 
relativistic in heavy atoms, the EDM of the atomic systems becomes larger than the 
electron EDM d e . In fact, the EDM of Cs atom has a large enhancement factor, namely 

d Cs ^ 914 [nunnnmsDEiEiiiBiiin]. 

However, it is also believed that the electron EDM might well be rather small 
compared to the neutron EDM. This is based on the EDM operator which is derived 
from the supersymmetry model calculations [21111^11221121111211 1221- The EDM interaction 
Tiedm is written for ipi fermion field with the EDM coupling constant di, 

n edm = —dSiV^^iF^ (1.2) 

where F^ u denotes the electromagnetic field strength. Up to now, there are many efforts 
to determine the strength of the EDM coupling constant di from the supersymmetric 
model calculations. 

At the present stage, however, it is not clear how large the electron EDM d e and 
neutron EDM d n should be. But most of the estimations of the EDM suggest that the 
electron EDM must be much smaller than the neutron EDM [TIJ. This means that it is 
better if one can measure the neutron EDM from atomic systems. 

Due to the presence of Schiff's theorem, it has been long believed that the EDM 
of nuclear systems should be quite small, and one obtains the nuclear EDM from Schiff 
moments jHj. 

However, the definition of the Schiff moment does not seem to be well established. 
Sushkov, Flambaum and Khriplovich jH] first introduced the Schiff moment which is 
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generated by P- and T- violating nucleon-nucleon interactions ^U]> and in this case, the 
Schiff moment is not directly proportional to the individual EDM of nucleons. However, 
according to the definition of the Schiff moment by Khriplovich and Lamoreaux [11J, 
the Schiff moment S is described in terms of the individual nucleon EDM d PiU as 

where Rq and I denote the radius and the total spin of the nucleus, respectively. 

This Schiff moment S can be obtained from the unitary transformation of the 
EDM Hamiltonian in terms of exp(i J2t=i • Pi) where Pj denotes the nucleon 
momentum operator in nucleus. But we should note that this type of the unitary 
transformation makes the calculation more complicated than the calculation with the 
original Hamiltonian. This is entirely due the fact that the nuclear system is not a 
neutral system, and therefore there is no Schiff shielding in the finite size effects of 
the atomic EDM. Therefore, even though the unitary transformation that translates 
the nucleon coordinates by the nucleon EDM d^r cancel out the EDM terms of the 
Hamiltonian that couple with the external electric field, it induces new EDM terms in 
the Hamiltonian some of which correspond to the Schiff moment. In this case, however, 
the translation of the nucleon coordinates affects the nuclear wave functions since the 
unitary transformation should naturally change the nuclear wave functions as well, and 
in fact, one can easily see that the EDM terms of the Hamiltonian that arise from this 
unitary transformation can be canceled out by the EDM terms from the nuclear wave 
function effects. Thus, the procedure of the unitary transformation should not change 
the EDM evaluation from the one that starts from the original EDM Hamiltonian. 
Instead, one should be careful for the evaluation of the nuclear EDM with the unitary 
transformation. If one solves the eigenvalues of the transformed Hamiltonian properly, 
then one finds the same result as the one from the original Hamiltonian. 

Also, we should add a comment on Schiff 's theorem why the unitary transformation 
can give a proper description of the EDM. In Schiff 's theorem, the unitary 
transformation cancels out the EDM terms completely from the Hamiltonian. In this 
case, the wave functions without the EDM terms can be employed since they are the 
eigenstates of the unperturbed Hamiltonian without the EDM terms. 

This point will be discussed in detail in ref. j2Hj, and here we only note that 
the EDM Hamiltonian of the present treatment is the same as the original Hamiltonian 
before the unitary transformation of the EDM operators. Further, we evaluate the EDM 
contributions corresponding to the EDM from the Schiff moment. But it turns out that 
this EDM contribution from the Schiff moment is very small as will be discussed in 
detail later. 

In this paper, we present a new calculation of the atomic EDM starting from the 
microscopic interactions in atomic and nuclear Hamiltonian. Here, we calculate the 
atomic EDM by considering all possible cases of the second order perturbation energies 
which are proportinal to the individual particles EDM dj and the external field E ext . 
Among many terms of the second order perturbation energies, the Schiff shielding can 
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be viewed as the cancellation between the first order perturbation energy of diE ext , and 
the second order perturbation energy of the two terms (djVjV^) and (er • E ext ) where V c 
denotes the Coulomb force of the corresponding system. 

The evaluations of the EDM are carried out in the three steps. In the first step, 
we take the nucleus as a point charge, and thus the finite size effects are neglected. 
In this case, we make the first and the second order perturbation theory and confirm 
that the first and the second order perturbation energies cancel out each other, and 
therefore there is indeed no EDM observable left due to Schiff's theorem. As is well 
known, this Schiff shielding can be obtained by the unitary transformation of the total 
atomic Hamiltonian for the point nucleus case. 

In the second step, we calculate the second order perturbation theory by taking 
into account the finite nuclear size effects. In this case, however, we consider atomic 
excited states as the intermediate states, but the nucleus is in the ground state due to 
the operators involved in the processes. From this calculation, we obtain the atomic 
EDM ((1a) which can be described by 

d A ^d n ~ -0.7 x l^Z 2 AU n (1.3a) 

for the nucleus with Z protons and A nucleons. Here Rq, Oq and AE e denote the 
nuclear radius, Bohr radius and the average electronic excitation energy, respectively. 
This result of eq.(1.3) should correspond to the part of the atomic EDM which arises 
from the Schiff moment, but the microscopic calculation shows that it is, in any case, 
extremely small. 

In the third step, we calculate the second order perturbation energy of EDM in 
nucleus, and we consider nuclear excited states as the intermediate states, but the 
electrons are in the ground state due to the operators involved in the processes. This 
is a new mechanism that generates the second order perturbation energy of the atomic 
EDM. Indeed, at a glance, one may feel that one should obtain a very small effect on the 
EDM from the nuclear excitations which are always of the order of 10 MeV. However, 
the operators are also written in terms of the nuclear variables and therefore they can 
give rise to a normal nuclear perturbation result. As one knows, the physical quantities 
in the perturbation theory in nucleus are usually not very small. In fact, we can make 
a rough estimation of the atomic EDM d A as 

Ze 2 

d A ~ a 77i p dn ~ 0.03Zd n (1.36) 

IXhiAtto 

where we made use of the following expressions for energy denominator AEa and nuclear 
radius R , 

41 

AE A ~ MeV, R ~ 1.2A 1/3 fm. 

From the realistic shell model calculations, we obtain the EDM of the deuterium, Xe 
and Hg atomic systems as 

d D ~ 0.017d n (1.4a) 



Microscopic Calculations of Nuclear EDM 



5 



d Xe ~ 1.6d n (1-46) 

d Rg ~ -2.8d n . (1.4c) 

These numbers are surprisingly large compared to eq.(1.3), and one sees that the nuclear 
EDM can be observed as the outer nucleon EDM. Also, the reason why one obtains 
rather large values for the heavy nucleus must be due to some coherence of the Coulomb 
interactions. Roughly, the EDM of the heavy nucleus is proportional to Z as shown in 
eq.(1.3b). 

From the above numbers together with the experimental observation of the atomic 
EDM G?Hg, we obtain the most severe constraint on the neutron EDM d n as 

d n ~ (0.37 ± 0.17 ± 0.14) x 10~ 28 e • cm (1.5) 

where we employed the observed values of the atomic EDM <i Hg 

d Ug ~ -(1.06 ± 0.49 ± 0.40) x 10~ 28 e ■ cm. (1.6) 

The EDM of neutron d n in eq.(1.5) is almost three orders of magnitude smaller than 
the direct measurement of eq.(l.l), and, even more, it seems to suggest that the d n is 
already quite close to a finite number. 

We believe that the finite number of the neutron EDM is quite reasonable. Since 
the size of neutron is ~ 10~ 13 cm and the parity violation should amount to ~ 10~ 6 , 
therefore, the pure T-violation is ~ 10~ 9 which is close to the baryon-photon ratio where 
the fundamental symmetries must have been violated. 

This paper is organized in the following way. In the next section, we treat the 
deuterium case, and clarify why the nuclear perturbation theory can give rise to the 
appreciably large EDM after Schiff's shielding effect is removed. In section 3, we present 
a general formalism of the EDM in nuclear systems with atomic electrons and show that 
the electric dipole operators with nuclear variables can contribute to the atomic EDM, 
and we apply the present formalism to the Xe and Hg atomic systems. Also, we show the 
reason why the atomic EDM due to Schiff moment is very small. Section 4 summarizes 
what we clarify in this work. 

2. EDM of deuterium 

In order to measure the EDM of atomic nuclei, we first have to face Schiff's theorem. 
The theorem states that the EDM operator of the particle may not have any effects on 
the physical observables since the effect of the EDM operator may be absorbed into the 
Coulomb interaction due to the translational property of the interaction Hamiltonian. 
In this case, one cannot measure any effect due to the EDM operator, and this is called 
Schiff shielding. 

However, the finite size of the nucleus brings about the EDM operator effects on 
the spectrum. In order to see the finite size effects in nucleus, we have to carefully write 
the Hamiltonian such that the effects can be seen explicitly. 
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Here, we first write the Hamiltonian for the deuterium system with explicit nuclear 
variables. The unperturbed Hamiltonian H of the deuterium system can be written 

where m and M denote the masses of electron and nucleon, respectively, p and r are the 
momentum and coordinate variables of electron while R p and R n denote the proton and 
neutron coordinates, respectively. R and P are the relative coordinate and momentum 
of the proton and neutron and are defined as 

R Rp Rn P 2^"^*^ P«-)' 
We choose the origin at the center of the deuteron, and thus we set 

Rg = + = ^' 

Also, Vjvjv(|R|) denotes the nucleon- nucleon interaction. 

Now, the perturbative Hamiltonian due to EDM can be written 

ed e -(r-R p ) ed p ■ (r - R p ) ed n • (r - RJ 

-Hp," 



'■edm 



|r-R p | 3 |r-R p | 3 |r-R n | 3 



+ ed ^ R - (d e + d p + d n ) • E ext + e(r - R p ) • E ext (2.2) 

where d e , d p and d n denote the intrinsic EDM of the electron, proton and neutron, 
respectively. They are related to the spin operator as 

dj = di<r i 

which can be obtained by the nonrelativistic reduction of eq.(1.2). The last term of 
Hedm comes from the normal electric dipole moment of the whole deuterium system 
which couples with the external electric field E ext . 

Eq.(2.2) can be rewritten in terms of the variable R 



Hp dm. 



ed e • (r - ±R) ed p • (r - ±R) ed„ • (r + ±R) 



edm '" |r - IR|3 " " |r - |R|3 | r + ±R| 3 

+ ed^_R _ ^ + ^ + d ^ _ ^ + e ^ _ 1 R ^. E ^ (2 3) 

2.1. Finite size of nucleus 

The atomic orbit is quite far from the nuclear radius, and therefore we can expand in 
terms of y. The unperturbed Hamiltonian becomes 

H ° = ^-j + w +VNN{m (2 - 4a) 
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where S in eq.(2.4b) is defined as 

S= 5 --^cos 2 Q (2.4c) 

where 6 denotes the angle between the electron coordinate r and the nucleon coordinate 
R. The perturbative Hamiltonian H e d m can be written as the term arising from the point 
charge H^ m together with the finite size Hamiltonian H™ np to the order of (R/rf 

= - e(de + d ; 3 +dn) ' r - (de + d p + d n ) • E ext + er ■ E ext (2.5a) 

H ^ = + d P + d«) • r + - ^eR • E ext (2.56) 

2.2. Schiff's shielding 

Now, we can calculate the EDM of the deuterium system from all possible cases of the 
perturbation energies up to the first order in dj and E ext . We write the ground state 
wave function of the deuterium system as 

^ = V4 O) (R)®0i O) (r) (2.6) 

where -?/>d(R) and 4> e (r) denote the wave function of deuteron and electron, respectively. 

First, we evaluate the first order perturbation theory and obtain the EDM energy 
of the deuterium AE^ 

AE^ = -(d e + d p + d n )E ext . (2.7) 
The second order perturbation theory without the finite size effect can be written as 



^gV(°)|(d e + d p + d n ) ■ V^o(r)|^ ) ^))(^ ) 0("-)|e^ |^g ) 

E nD ,n e — Eq 



E 



(<0( o )|e^|^V("-))(^ ) 0("-)|(d e + d p + dn) ■ VA (r)|< 

rp rp 

(2.8) 

g 

where Aq(y) = - is introduced and E denotes the energy eigenvalue of the ground state 
r 

in the deuterium system. Here, we should note that the nuclear part is always in the 
ground state and has no effect. 

Now, we make use of the following identity 

VA)(r) = tpA (r) = i[p, A (r)} = --[p, H ] (2.9a) 

e 

Hol^nl) = Eol^l), H \¥£) = E nD , ne \¥£). (2.9b) 
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Therefore, we can sum up all the intermediate states since the energy denominator 
En D ,n e ~ E Q cancels out, and thus obtain 

AE% = (d e + d p + d n )E ext . (2.10) 

This is exactly the same as the first order perturbation result with opposite sign, and 
thus it cancels out the first order perturbation result of the EDM in deuterium. This is 
just Schiff's shielding. 

2.3. Finite size effect (atomic excitation) 

Now, we consider the second order perturbation theory with the finite size effects of the 
nucleus taken into account. Here, we first calculate the atomic excitation in the same 
way as the previous subsection, but consider the first term in eq.(2.5b) together with the 
last term of eq.(2.5a). The second order EDM energy with the finite size effect becomes 

AEf] = -2 £ 1 



x <<#l^(d e + d p + d n ) • r|^V^)(^ D Vi ne) |e^ e:ct |^ Vi 0) ) (2.H) 

where the intermediate states I^d^^J 1 ^} denote the nuclear and atomic excited states. 
From the operator properties, it is clear that the nuclear state should be in the ground 
state in the intermediate state. 

We can evaluate eq.(2.11) by making use of the closure approximation, and here 
we take < E QtrLe — E >~ 10 eV as an optimistic value. In this case, we obtain 



AE ?s = F : ^r~(de + d p + d n )E ext 



e 



2 



x(^|^ 2 |^)(0i°)|^|^)). (2.12) 



, cos 2 6 



We can calculate eq.(2.12), and obtain the second order EDM energy and the EDM of 
the nucleus 

AEf] ~ -1.6 < S > xl0~ 9 (4 + d p + d n )E ext (2.13a) 

d D ~ -1.6 < S > xl0- 9 (d e + d p + d n ) (2.136) 

where < S > denotes the expectation value of the S, and if we calculate < S > with 
the density distributions with the spherical symmetry, then < S >= 0. Therefore, the 
term which corresponds to the Schiff moment vanishes. Since < S > is zero or of the 
order of unity, this dp is extremely small, and there is little chance to observe it. 
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2.4- Finite size effect (nuclear excitation) 

From eqs.(2.5), one sees that the second order perturbation energy of EDM has a 
contribution from the nuclear excitations. In this case, the atomic state is kept in 
the ground state, and we obtain the second order EDM energy as 

n D ,n e En D ,n e ~~ -^0 

where the intermediate electron state must be in the ground state due to the orthgonality 
condition between the ground state and the excited state since the operators do not 
depend on the electron coordinates. 

The excitations of deuteron are from continuum states. In this case, we can describe 
the state \tia) and the excitation energy E nA as 

|Vi nD) ) = ^e X p(ik.R) (2.15a) 

E nD = ^- (2.156) 

Therefore, we can write eq.(2.14) as 

AEfJ = e 2 d n E ext J ^dWR'^^M^MR')^ 

x cos 6 cos 6' exp [ik ■ (R - R')] (2.16) 

where the spin part of the deuteron wave function is already carried out. Also, we took 
only the S-state of the deuteron wave function, k B is related to the binding energy of 
the deuteron, and can be written 



k B = ^ME B ~ 46 MeV/c. (2.17) 
The integration over the momentum k can be easily done and we obtain 

AEf] = e 2 d n E ext ^- J d 3 Rd 3 R' eXp [ ~^|* ~ ^h oWM*)^ cos Acqs ff. (2.18) 
The angular integrations can be also done, and we obtain 

AEfJ = ^-e 2 d n E ext Mk B J R 3 dRdR'^ D (R)^ D (R') 
( cosh(k B R<) smh(k B R < )\ ( 1 1 \ , , „ N . « 

x ~ J + <w) e * pM) (219) 

where R > (R<) denotes the larger (smaller) one of R,R'. 

For the numerical evaluation, we use a simple deuteron wave function of the 
following shape 

MR) = (^) 2 exp(-/3i?) (2.20) 
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where we took the value of (3 as (3 = = 46 MeV/c. Finally, we obtain the second 
order EDM energy and the EDM of the nucleus 

AE$ ~ 0.0l7d n E ext (2.21a) 

d D ~ 0.017d n . (2.216) 

This is not at all a small number, and we believe that the EDM measurement of the 
deuterium system can be well competing with the neutron EDM. 

In fact, there is a proposal of EDM collaboration to measure the EDM of deuterium 
|27j . According to their proposal, they can measure the deuterium EDM up to 

do ~ 1CT 27 e • cm. 

In this case, the neutron EDM d n becomes 

i~5x 1(T 26 e • cm 

which is just competing with the present upper limit of the neutron measurement of 
eq.(l.l). 

At this point, we should discuss the relation between the EDM [eq.(2.14)] due to the 
nuclear excitation and the EDM arising from the T-, P-odd nucleon-nucleon interactions. 
The T-, P-odd nucleon-nucleon interactions should be obtained from the second order 
Feynman diagram of the color EDM of the quark-quark-gluon vertex operator together 
with the normal quark gluon vertex terms. Therefore, it is not clear how one can relate 
this T-, P-odd nucleon-nucleon interaction to the neutron EDM. Further, this T-, P- 
odd nucleon-nucleon interactions should be quite complicated since they are two-body 
nuclear interactions. In this case, the calculation of the EDM should be much more 
involved than the present method, and the EDM arising from the T-, P-odd nucleon- 
nucleon interaction is an indirect evaluation of the atomic EDM compared to the present 
calculation. The EDM due to T-, P-odd nuclear forces is discussed by Flambaum, 
Khriplovich and Sushkov [2H], and they obtained the EDM of the deuteron as 

d D = 0.2 x 10~ 21 £ e • cm 

where £ is a parameter related to the T-, P-odd nuclear forces. We believe that the two 
different methods should give the different contributions to the deuteron EDM. 

3. EDM of Xe and Hg atomic systems 

Up to now, we have discussed the deuterium which is the simplest nucleus with one 
electron. Now, we discuss more complicated nucleus, and here we treat the Xe and Hg 
atomic systems since there are some measurements of the EDM in these atomic systems 
jH EH El and also there is a proposal to measure the EDM of the atomic and nuclear 
system [29J. 
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In the same way as the deuterium case, we first write the Hamiltonian of the total 
atomic and nuclear system. 

The unperturbed Hamiltonian H of the Xe system can be written 



#o = E 



Pi ^ 



2m 



1 \ r i Rj 



1 Z 



r, - r,- 



+ E^ + 2E^(|R l -R,| 



1 



E 



(3.1) 



^ 2M _ 

where r^, denote the coordinate and the momentum of the electron while Rj, Pj 
denote the nuclear variable and momentum, respectively. 

On the other hand, the perturbed Hamiltonian coming from the EDM is written as 



z z ed* ■ (r 

n edm / j / j 

i=l j=l 



r,-R,f 



Rj) , VV ed e-^ 



i=i j ^i 



i=ij=i I 1 ' 
z 



A z edV • (Ri - Rj) 



EE 

i=i ^ 



R: 



(3.2) 



- E d e • E ext - E d JV • Eezt + e E( r * " R») " E e 

i=l i=l i=l 

where the summation over Z in nucleus means that it should be taken over protons. 
The EDM of the nucleon can be expressed in terms of the nucleon isospin as 



n 2 



;i+t^ + (i-^' 



3.1. Finite size of nucleus 

Now, we evaluate the finite size effects on the second order EDM energy in heavy nucleus. 
The unperturbed Hamiltonian of the atomic and nuclear system becomes 



#o = E 



i=\ 



2m V; 



r, - i\ 



R; — R, 



Here, we ignore the finite size effect of the unperturbed Hamiltonian. 



Now, the perturbed Hamiltonian H^ m from the point charge and the Hamiltonian 



(3.3) 



Hfdm w ^h the finite size can be written up to the order of (Rj/ri)' 



H 



(0) 
edm 



-E 



eZd l e • -1 

rf 



A 

Eed 

u=i 



_ edj • (rj 

3 |„. 



r 



.13 
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/ Z A \ Z 

~ E d e + E d V -E^ + e^ivE^ (3.4a) 

\i=l i=l / i=l 



ttUs) = V- 

edm / y 



Z 



Z A 

2 
3 

3=1 3=1 



i=i 

1 Z fRi-R,-) e 1 



e 

^5 



"EE • - 9 E(! + • E««* (3-46) 

i=i 3 yi l 1 ^ ^Vj'I z i=i 

where Sji in eq.(3.4b) is defined as 

5 15 

Sn = 2 - y cos2 e ^ ( 3 - 4c ) 

where 6jj denotes the angle beween the electron coordinate and the nucleon 
coordinate Rj, and can be given as 

cos Qji = sin 9j sin Q { cos(0j — (pi) + cos 9j cos 0j. 

3.2. Schiff shielding 

When we treat the nucleus as a point particle, the first order and the second order EDM 
energies cancel out each other. That is, we first estimate the first order perturbation 
energy of the second term — (j2f=i d* + J2f=i d^r) • E ext in eq.(3.4a), and indeed this is 
finite. Then, we evaluate the second order perturbation energy of the first term and the 
last term in eq.(3.4a). This second order perturbation energy exactly cancels out the 
first order perturbation energy, and this is indeed the Schiff schielding. The calculation 
in detail is just the same as the deuterium case, and therefore we do not repeat the 
procedure of Schiff shielding. 

3.3. Finite size effect (atomic excitation) 

Now, we consider the second order perturbation energy with the finite size effects of the 
nucleus. In order to evaluate the second order EDM energy, we first need to have the 
atomic and nuclear wave functions. We write the wave function of the total atomic and 
nuclear system by 

*t'" e) = V4 nA) (Ri, • • • , Ra) ® 0i ne) (ri, • • • , v z ) = ^ A) € £) (3-5) 

where ua, n e denote the quantum number of the excited states, and n\ — 0, n e = 
mean the ground states of the systems. Here, we assume that the atomic state has the 
ground state with spin zero while the nuclear ground state has one outer neutron with 
spin \. This is mainly because we consider the 129 Xe and 199 Hg atomic systems in this 
paper. 

The second order EDM energy with the finite size effect becomes 

AE fs = - E p 2eE T F (^V^lgigl^^)^^!!:^!^^) 

n A ,n e ^n A ,n e j = l 

(3.6) 
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where E denotes the ground state energy of the whole system. From this equation, 
it is clear that the intermediate nuclear state is the ground state = due to the 
orthogonal condition since the operators only involve the electron coordinates in the 
last term of eq.(3.5). Further, H^2,'° is defined as 



edm / j 



i=l 



Z A 

I 2 



(I! • r, £ S 5i R) I v, • X & N SjiR) 



^5 



(3.7) 



Eq.(3.6) can be evaluated by employing the closure approximation 

AE?}^ 2&2Eext x 




< E 0t n e — E > 
Z A 



d* • r*) J2 SjiR* + E( d ^ ■ **) V§ ] 3|0>. (3.J 

i=i i=i 



This becomes 



AEfJ ~ - - ^^ e g , ^L^l^i^li (3.9a) 



l{2) _ 2e 2 ZE ea;f rf n < % cos 2 8j > (R 2 ) 
< Eo,n e — E > 

where the first term in eq.(3.8) vanishes since the spin of the atomic states is assumed 
to be zero. 

Here, a denotes the Bohr radius of the atomic system and can be written as 
a = „ 2 . Therefore, the atomic EDM from the atomic excitations becomes 



d 



A 



2e 2 Z d n < SjiCcxPOi > (R 2 ) 



mZ 2 e i <2q 

~ -1.0 x 1(T 9 < S jt cos 2 Bi > Z 2 Ald n (3.96) 

where we take < E Q ^ — E >~ m(Ze 2 ) 2 and (R 2 ) ~ r\Ah with r = 1.2 fm. 
< Sji cos 2 6i > is zero if we evaluate it with the density distribution with the spherical 
symmetry. Since < Sji cos 2 &i > is zero or of the order of unity, this d& is very small 
and there is again little chance to observe it. 

Eq.(3.9) just corresponds to the EDM of the atomic system calculated by the Schiff 
moment in [TT] even though there must be some difference by a factor of 2. 



3.4- Finite size effect (nuclear excitation) 

Now, we consider the second order EDM energy due to the intermediate nuclear 
excitations. The important point is that we should estimate all possible cases of the 
second order processes with the condition that the EDM dj of the individual particles 
and the external field E ext must be in the first order. This process arises from the finite 
nuclear size effects in the EDM Hamiltonian. The second order EDM energy can be 
written as 
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X (^TVH E ^ [(1 + Tf )^ + (1 - Tf )^] • (1+ |^; 3 Rj) |^ 0) ^ 0) ) (3-10) 



|Ri-R,f 



From eq.(3.10), it is clear that the electron state should be in the ground state (n e = 0) 
due to the orthogonality condition since the the operators in the first term in eq.(3.10) 

A 

have only the nuclear coordinates. Here, the relation E R» = should be taken since 

we set the center mass coordinate to zero. 

Here, we carry out a rough estimation of eq.(3.10) before the realistic shell model 
calculations. Eq.(3.10) may be written as 

AEf] -^<R> E ext d n ^^ (3.11a) 

Therefore, the EDM oIa of the total atomic system can be written as 

Ze 2 

where we made use of the relations AE ~ huo ~ MeV, R ~ 1.2A 1 / 3 fm. 
3.5. Xe atomic system 

Now, we calculate the 129 Xe case in which we assume a simple single particle shell model 
state, namely, 3si for neutrons together with the 2d,3 <g> 2 + . Further, the atomic states 
stay in the ground state, and therefore the electron wave functions are not written here. 
Thus, we write the nuclear wave function as 

\Va*) = Vl-a 2 \u(3si) : h + a\v(2ds) ® 2+ : J) (3.12a) 

2 Z 2 2 

where a is a parameter which should be determined from the realistic calculations. 

In this case, the intermediate states that contribute to the second order EDM energy 
(eq.(3.9)) for f(3si) are restricted to the following two states 

|n) = \u(3pi) : \) , k[(3 S i) 2 [0 + ], (2p^} : \) (3.126) 

where the second term comes from the hole state contribution. Here, we note that only 
the above two states can contribute to the matrix element since the matrix elements of 
the radial component R from higher states vanish. 

Further, we assume that the 2 + state does not contribute to the atomic EDM since 
it is a collective state and the operators we are considering should not give a large 
contribution to the second order perturbation energy. This is mainly due to the fact 
that the two operators have different rank tensors which generate 1~ and 0~. 

Therefore, we take the following intermediate states corresponding to the second 
term of eq.(3.12a) 

\n) = \u(3 Pl ) ® 2 + : i) , |i/[(2d f ) 2 [0 + ], (2p f H 2+:^). (3.12c) 
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With these simplified configurations, we can calculate eq.(3.12) and find the second 
order EDM energy due to the nuclear excitations, 



a 7-1(2) _ Ze 2 d n E ext 



UJ 



;i-a 2 ){M3 S p|Z|K3 P| ))M3^)|^IK3 S i)) 



-M3si)|^IK2Pi)>M2Pi)I^IK3si)>} 

{1 11 U Wj 1 

(u(2ds 2 ) ® 2+ : -|Z|zy(3p f ) ® 2 + : -> (i/(3p f ) ® 2+ : -| — |„(2d f ) ® 2 + : -> 

- (z/(2da) ® 2+ : ^|^(2p f ) ® 2 + : i) (i/(2ps) ® 2+ : I|^V(2d § ) ® 2+ : i) 

(3:i3) 

where a is taken to be a 2 ~ 0.5. Further, in evaluating eq.(3.10), we made an 
approximation of |Rj — Rj| 3 i? 3 which may induce some error of 20 ~ 30 percents in 
eq.(3.13). 

Here, we take the harmonic oscillator wave functions for the single particle state, 
and uj is given as uj ~ MeV. The matrix elements can be calculated as 

A3 



M3 S ,)|Z|,(3^)) = i^ (3.14a) 
M3^)|^IK3 S| )) = ^M. (3.14o) 



(i/(2d f ) ® 2 + : ±|Z|i/(3p f ) ® 2+ : i> = 1^ J- (3.14c) 
(i/(3p f ) ® 2+ : i|^|K2rf f ) ® 2+ : I) = ^sl^Muj ^.Ud) 



^)\ZW7p$) = -\^ (3.14e) 

WtW 3 ^ (3 - 14/) 



M2d.) 2+ : 1|Z|K2P|) ® 2+ : I> = -^/^- (3.14,) 

M2p f ) ® 2 + : I|^>(2ds) ® 2+ : I) = g |/|m,. (3.14/,) 

Therefore, we can evaluate eq.(3.13) and obtain dx e 

d Xe ^ l-6d n . (3.15) 

This is a surprisingly large number, but we believe that the number must be reliable 
within a factor of two or so. 
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For the proton EDM, we should consider the core excitation, and this is rather 
difficult to estimate since the simple-minded calculation cannot give a finite number. 
This is simply due to the operators that induce the intermediate states. If we start 
from the core which has the + state, then the dipole operator r?Rj • E ext can only 
generate the isovector 1~ states while the interaction part can only generate 0~ 
states. Therefore, if we start from the core with + state, then there is no contribution 
to the atomic EDM. 



3.6. Hg atomic system 

Next, we calculate the 199 Hg case in which we also assume a simple single particle shell 
model states as 

1 1 X 

\*Ae) = \A - «i - «2 K3pi) : -)+a 1 \v(2f,)®2 + : -) + a 2 \v(?>pz) ®2 + : -> (3.16a) 

2 Z 2 Z 2 z 

where a.\ and a 2 are taken to be a\ ~ 0.25 and a\ ~ 0.25. Here, we consider the 
following intermediate states for the first term of eq.(3.16a), 

|n) = |i/(4*i) : \) , |i/[(3pi) 2 [0 + ], (3s^] : \) (3.166) 

where the second term comes from the hole state contribution. 
For the second and the third terms of eq.(3.16a), we take 

|n) = |i/(3dj) ® 2 + : I) , |i/[(2/.) 2 [0 + ], (2d § H ® 2 + : I) (3.16c) 
|n) = |i/(3ds) <g> 2+ : i) , |i/[(3p s ) 2 [0+], {2d*)- 1 } ® 2 + : -). (3.16d) 

2 2 2 2 2 

(2) 

In this case, A_Ej s can be written as 



^gi(2) _ Ze 2 d n E ext 



[l-a\- al) {M3pi)|^IK4si)>M4si)|^IK3pi)> 



-M3pi)|Z|K3 S i))(K3 S i)|^IK3pi))| 
+« 2 | (i/(2/.) ® 2+ : i|Z|i/(3d.) ® 2 + : \)(v{*d*) ® 2+ : i|^V(2/ f ) ® 2+ : i) 
- M2/ f ) ® 2+ : \\ZH2ds ) ® 2+ : ±)(i/(2d f ) ® 2+ : I|^?>(2/ s ) ® 2+ : ±>j 
+a 2 |(K3P3) ® 2+ : ^|Z|i/(3d|) ® 2+ : ^>(i/(3ds) 2+ : ±|^>(3ps) ® 2+ : ±) 
- (i/(3p f ) ® 2+ : i|Z|i/(2d f ) 2+ : i>(i/(2d f ) ® 2+ : I|^2>(3 P| ) ® 2+ : i> 



(3.17) 
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Here, the matrix elements become 



^)\ZH^)) = -\^ (3.18a) 
Ws k )\^\^ n) ) = ^=M. (3.186) 



(1/(2/3) ® 2+ : i |Z|i/(3d f ) ® 2+ : i) = ^ (3.18c) 

(i/(3d f ) ® 2+ : ^l^?IK2/|) ® 2 + : ~> = || J| M, (3.18d) 



(i/(3p 3 ) (8) 2+ : -|Z|i/(3d 3 ) ® 2+ : -) = —J-^— (3.18e) 

X u 37 

(i/(3d3) ® 2+ : -|— — — |zy(3p3) ® 2+ : -) = =Mw (3.18/) 



M3^)|^IK3^)) = ^^Mc (3.18/,) 



(i/(2/j) (8 2 + : ^|K2rf|) ® 2+ :\) = ± ^ -L_ (3.18i) 
(K2d|) ® 2+ : ^l^k(2/|) ® 2+ : i> = A j| Mu; (3.18/) 



(z/(3p.) ® 2+ : i|Z|i/(2d.) ® 2+ : ~> = (3.18k) 

(v(2ds) ® 2 + : I|^|„(3pj) ® 2+ : I) = ^Ma;. (3.18/) 
Therefore, we can evaluate eq.(3.15) and obtain dn g 

d Hg ~-2.8d n . (3.19) 

3. 7. Neutron EDM d n 

Since we have now the relation between the atomic EDM and the neutron EDM d n , we 
can deduce the constraint on the d n from the experimental values of dx e and dug- In 
fact, they are given as |U H] 

d Xe ^ (0.5 ± 2.4 ± 0.1) x 1(T 27 e • cm (3.20a) 



d Hg ^ -(1.06 ± 0.49 ± 0.40) x 10~ 28 e • cm. 



(3.206) 
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By making use of eqs.(3.15) and (3.19), we obtain the constraints on the d. 



71 



d n ~ (0.3 ± 1.5 ± 0.1) x 10~ 27 e • cm 



(3.21a) 



d n ~ (0.37 ± 0.17 ± 0.14) x 10~ 28 e • cm. 



(3.216) 



This is quite a stringent constraint on the neutron EDM d n even if we compare it to 
the direct measurement 1^1 E3 of the neutron EDM d n in eq.(l.l). In particular, the 
neutron EDM d n from the Hg measurement seems to suggest that the d n is getting quite 
close to a finite number. It should be very interesting if one could measure the EDM of 
Xe nucleus up to the same order as the Hg case. We can expect a finite neutron EDM. 

Since the neutron EDM d n is smaller than the value given in eq.(l.l), this may 
give rise to some serious problems for the supersymmetry model calculations [21] which 
already predict the neutron EDM d n comparable to the value of eq.(l.l). There must 
be some suppression mechanism |3Uj . 

4. Conclusions 

We have carried out the microscopic calculations of the second order perturbation energy 
of the atomic EDM for deuteron, Xe and Hg atomic systems. The basic point of our 
calculation is the new mechanism that generates the atomic EDM where the atomic 
states stay in the ground state while the nuclear excitations arising from the nuclear 
dipole operators are taken into account. All the calculations are done with the shell 
model wave functions even though we have employed a simple-minded shell model state. 
This may induce some errors of a factor of two or so, but it is for sure that the EDM of 
Xenon nucleus dx e is comparable to the neutron EDM d n itself. 

The atomic EDM that arises from the atomic excitations is related to the Schiff 
moment where the atomic EDM is proportional to the mean square radius of nucleons, 
< R 2 >. It turns out, however, that the EDM due to the Schiff moment is very small if 
we evaluate it microscopically. Therefore, the atomic EDM from the atomic excitations 
is practically impossible to measure. 

From the present work, one sees that the atomic EDM should be appropriate for 
the experimental observation of the neutron EDM. Even for deuterium or D 2 molecule 
which is better since one gets a factor of two compared to the deuterium, the EDM of 
deuterium is a few percent of the neutron EDM. We believe that this number can well 
compete with the neutron EDM measurement which is usually quite difficult due to the 
short life time of neutron. 
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